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u.  abstract  (M*Kmvm2oo womsf  ihe  effects  of  dynamic  stall  lift 
stability  of  a  hingeless  rotor  are  investigated,  both  experimentally  and  analytically 
The  eephasis  is  cn  the  oorrelaticn  with  measured  regressing  lead-lag  mode  darping 
levels  of  a  soft-inplane,  three-bladed  model  rotor,  coerated  untrimed.  The  correlate 
covers  awide  range  of  test  conditions  for  several  values  of  rotor  speed,  collective 
pitch  angle,  shaft  tilt  angle  and  advanve  ratio.  It  includes  cases  that  vary  from  nea 
zero  thrust  conditions  in  hover  to  highly  stalled  forward  flight  conditions  with 
advance  ratios  as  high  as  0.55  and  shaft  tilt  angles  as  high  as  20  .  Both  the  expep.- 
rental  and  analytical  blade  models  represent  a  sirple  model  of  a  hingeless  rotor  with 
rigid  blades  and  spring  restrained  flap- lag  hinges.  The  aerodynamic  representation  is 
based  cn  tie  CNERA  dynamic  stall  models  cerprising  virtually  independent  unified  lift 
ani  drag  node  Is.  He  nonlinear  equations  of  blade  nation  and  stall  dynamics  are 
perturbed  about  a  periodic  forced  response,  and  tie  damping  is  evaluated  by  the 
Floquet  eigenanalvsis.  In  ccnpariscn  to  the  linear  and  quasi  steady  stall  aerodynamic 
theories,  the  theory  with  dynamic  stall  lift  and  quasisteady  stall  drag  qualitatively 
improves  the  oorrelaticn,  and  adding  dynamic  stall  drag  provides  further  ouantitative 
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Summary 


The  effects  of  dynamic  stall  lift  and  drag  on  the  flap-lag  stability  of  a  hingeless  rotor 
are  investigated,  both  experimentally  and  analytically.  The  emphasis  is  on  the  correlation 
with  measured  regressing  lead-lag  mode  damping  levels  of  a  soft-mplane,  three-bladed 
model  rotor,  operated  untrimmed.  The  correlation  covers  a  wide  range  of  test  conditions  for 
several  values  of  rotor  speed,  collective  pitch  angle,  shaft  tilt  angle  and  advance  ratio  It 
includes  cases  that  vary  from  near  zero  thrust  condition  in  hover  to  highly  stalled  forward- 
flight  conditions  with  advance  ratios  as  high  as  0  55  and  shaft  tilt  angles  as  high  as  20° 
Both  the  experimental  and  analytical  blade  models  represent  a  simple  model  of  a  hingeless 
rotor  with  rigid  blades  and  spring  restrained  flap-lag  hinges.  The  aerodynamic 
representation  is  based  on  the  ONERA  dynamic  stall  models  comprising  virtually 
independent  unified  lift  and  drag  models.  The  nonlinear  equations  cf  blade  motion  and  stall 
dynamics  are  perturbed  about  a  periodic  forced  response,  and  the  damping  is  evaluated  by 
the  Floquet  eigenanalvsis.  In  comparison  to  the  linear  and  quasisteady  stall  aerodynamic 
theories,  the  theory  with  dynamic  stall  lift  and  quasisteady  stall  drag  qualitatively 
improves  the  correlation,  and  adding  dynamic  stall  drag  provides  further  quantitative 
improvement. 


Notation 


Unless  otherwise  stated,  the  symbols  below  are  nondimensional 
a  linear  lift  curve  slope,  ’■ad'1 

ad  damping  factor  for  dynamic  stall  drag 

b  blade  semi-chord,  (1/R) 

c  blade  chord,  tn 

Cj,  Ci  airfoil  sectional  drag  and  lift  coefficients 
Cu$  quasisteadv  stall  drag  coefficient 

Cdo  airfoil  constant  profile  drag  coefficient 

Ci0  airfoil  lift  coefficient  at  zero  angle  of  attack 

CT  thrust  coefficient 

Cti,  C-i  linear  and  quasisteadv  stall  lift  coefficient 

Cy  nonrotating  flap  damping  coefficient 

nonrotating  lead-lag  damping  coefficient 

D  airfoil  drag  force  per  unit  length  in  the  free-stream  direction,  {l/p  b  Q'R:) 

eh  blade  hinge  offset,  (1/R) 

c  phase  shift  parameter  for  dynamic  stall  lift 

ed  phase  shift  parameter  for  dynamic  stall  drag 

F-,  F,  airloads  per  unit  length  m  Dap  and  lead-lag  directions 

k  b/r 

kr  reduced  frequency.  Fig  4 

L  airfoil  lift  force  per  unit  length  normal  to  the  free-stream.  (1/p  6 

L0  apparent  mass  lift  normal  to  the  chord  line,  (1/p  4 

X  number  of  blade  segments 

P  flap  stiffness  parameter 

r  radial  distance  from  rotor  centei,  (1/R) 

rd  dynamic  stall  drag  frequency  parameter 


rotor  radtus,  in 

apparent  mass  parameter,  r  .d'1 

resultant  air  velocity  at  a  blade  section,  (1/flR) 

local  air  velocities  components  in  the  rotor  hub  plane,  (1/flR) 

dynamic  stall  lilt  frequency  parameter 

lead-lag  stiffness  parameter 

blade  radial  distance  measured  from  the  hinge,  (1/R) 

local  blade-fixed  coordinates 

flap-lag  stiffness  parameter 

blade  airfoil  section  angle  of  attack,  rod 

camber  of  the  airfoil,  rad 

quasisteady  stall  angle,  rad 

shaft  tilt  angle,  rad 

blade  flapping  angle,  rad 

lock  number 

circulation  per  unit  length,  UCi 

Circulation-like  drag  per  unit  length,  UC; 

width  of  the  ith  aerodynamic  element 

pitch-rate  coefficient,  rad'- 

airfoil  rotation  rate  with  respect  to  the  airmass 

blade  lead-lag  angle,  rad 

time  delay  parameter 

total  inflow  ratio 

blade  pitch  angle,  rad 

collective  pitch  angle,  rod,  dec 

advance  ratio 

negative  of  lead-lag  mode  damping  exponent,  sec  1 


Ci 
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rotor  solidity 

blade  section  inflow  angie.  Fig.  2 

nondimensional  time  or  blade  azimuthal  position,  rad 

frequency  term  due  to  blade  hinge  offset 

rotating  lead-lag  frequency,  (1/fl) 

nonrotating  flap  and  lead-lag  frequencies,  (1/11) 

rotor  angular  speed,  rad/ sec 

ith  blade  segment 

unstalled  and  stalled  components 

x  and  y  components 

dQ/dti 

equilibrium  state 
perturbed  state 
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1.  Introduction 


Very  little  is  known  about  rotor  blade  stability  under  dynamically  stalled  conditions 
in  forward  flight.  This  is  particularly  true  of  lead-lag  or  lag  modes  that  are  at  best  weakly 
damped  and  for  which  current  methods  of  predicting  damping  levels  merit  considerable 
refinement.  Given  the  complexity  of  the  flow  field  of  a  rotor  blade  in  dynamic  stall,  it  is 
imoortant  that  analytical  results  are  compared  with  a  broad  range  of  experimental  data. 
Accordmgly,  the  present  study  investigates  the  effects  of  dynamic  stall  lift  and  dynamic 
stall  drag  on  the  Cap-lag  stability  of  an  isolated  hingeless  rotor.  The  emphasis  is  on  the 
correlation  of  analytical  lead-lag  damping  results  with  experimental  data  from  tests 
performed  on  a  three-bladed  isolated  model  rotor  in  the  U.S.  Army  Aeroflighldynamics 
Directorate’s  7-  by  10-foot  (2.1  m  by  3.0  m)  wind  tunnel  (Ref.l).  The  experimental  model 
represents  a  simple  model  of  a  hingeless  rotor  with  rigid  blades  having  spring  restrained 
Cap  and  lag  hinges.  This  intentional  simplicity  helps  to  focus  the  testing  on  the  main  topics 
of  interest,  which  are  the  forward  flight  and  dynamic  stall  aspects  of  the  stability  problem 
The  structural  details  of  the  analytical  model  were  chosen  to  correspond  as  closely  as 
possible  to  the  experimental  model  to  remove  possible  sources  of  error  in  the  correlation 
The  database  comprises  lag  regressing  mode  damping  and  covers  a  wide  range  of  test 
conditions  for  various  values  of  rotor  speed  (£1),  collective  pitch  (0„),  shaft  angle  (as),  and 
advance  ratio  (*.)  It  includes  cases  that  vary  from  near  zero  thrust  conditions  in  hover  to 
highly  stalled  forward  flight  conditions  with  advance  ratios  up  to  0.55  and  shaft  tilt  angles 
to  20°. 

The  present  work  originated  as  part  of  a  research  program  to  validate  new 
analytical  methods  based  on  Floqaet  theory  for  predicting  the  effects  of  forward  flight  on 
aeroelaslic  stability  of  rotorcraft.  The  work  of  Ref.  2  on  flap-lag  dynamics  of  an  isolated 
blade  is  a  typical  early  application  of  Floquet  theory  to  solve  the  periodic  coefficient 
equations  that  arise  in  forward  flight.  Virtually  no  systematic  correlations  exist  that 
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address  the  validity  of  such  Floquet  theory,  even  though  it  is  now  m  relatively  widespread 
use.  The  Aeroflightdynamics  Directorate  undertook  a  carefully  designed  small-scale  model 
testing  program  to  obtain  data  for  such  validations.  The  results  of  the  experimental  tests 
are  reported  in  Ref.  1,  where  the  effects  of  periodic  excitations  m  forward  flight  presented 
no  particular  pr<  em  in  extracting  lag  damping  from  transient  response  data  Comparisons 
with  theory  are  reported  in  Refs.  3  4.  In  Ref.  3,  the  prediction  is  based  on  linear 
quasisteady  aerodynamic  theory  with  a  dynamic  inflow  unsteady  wake  theory  The  results 
were  satisfactory  only  for  very  low  values  of  collective  pitch  and  advance  ratio  and  it  was 
suspected  that  the  effects  of  airfoil  stall  might  be  responsible.  In  an  attempt  to  explain  the 
reasons  for  the  poor  correlation,  the  emphasis  of  the  investigation  shifted  to  the 
examination  of  the  influence  of  airfoil  stall  on  forward  flight  aecoelastic  stability  Earlier 
studies  (e  g.  Ref  5),  comparing  quasisteady  aerodynamic  theory  with  experimental  data  for 
the  hover  flight  condition  had  shown  that  airfoil  stall  strongly  influenced  flap-lag  stability 
in  hover  at  high  collective  pitch  angles  and,  furthermore,  that  quasisteady  nonlinear  airfoil 
stall  characteristics  could  successfully  account  for  these  influences  Consequently,  these 
effects  were  included  in  the  forward  flight  flap-lag  studies  in  Ref  4.  For  hover  and  very  iow 
advance  ratios,  the  anticipated  improvements  were  largely  realized,  however,  the  results  for 
high  advance  ratios  were  poor  and,  paradoxically,  in  some  cases  even  worse  than  the  linear 
theory  correlation.  As  a  result,  the  present  investigation  was  undertaken  and  a  nonlinear 
dynamic  stall  theory  was  included  in  the  analytical  predictions 

Rotor  blades  in  forward  flight  encounter  a  highly  complex  three-dimensional 
unsteady  flow  field,  which  typically  includes  regions  that  for  many  operating  conditions  are 
partly  stalled.  Several  dynamic  stall  models  have  been  developed  and  used  to  predict  time 
histones  as  well  as  damping  levels  (Refs.  6-10).  In  this  paper,  the  aerodynamic  description 
is  based  on  the  ON'ERA  dynamic  stall  unified  lift  and  drag  models  (Refs.  6  and  T)  The 
nonlinear  equations  of  blade  and  stall  dynamic:  Ci^  perturbed  about  a  periodic  forced 
response  to  compute  damping  levels  from  the  Floquet  eigenanalysis  The  emphasis  is  on 
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demonstrating  the  need  ter  including  dynamic  stall  in  the  prediction  of  lag  damping,  not  on 
the  relative  merits  of  different  dynamic  stall  models  or  of  the  Floquet  stability  analysis  vis- 
a-vis  other  stability  analysis  methods.  Given  the  breadth  of  the  study,  this  correlation  will 
provide  a  reference  comparison  with  other  correlation  studies  based  on  different  approaches 
to  modeling  dynamic  stall  and  to  predicting  stability  margins.  It  will  also  help  clarify  the 
effects  of  dynamic  stall  on  the  lag  stability  of  isolated  rotors. 

2.  Experiment  and  Data 

A  bnef  account  of  the  experimental  model  and  database  is  given  here  (for  details  see 
Ref.  1).  To  enable  the  use  of  a  simple  rigid  flap-lag  analysis  for  correlation,  the  three- 
bladed  rotor  used  flexures  to  simulate  hinged  blades  with  spring  restraint  and  coincident 
flap  and  lag  lunges,  the  effective  hinge  offset  was  0.11R.  The  untapered  rotor  blades  had  no 
pretwist  and  were  stiff  relative  to  the  flap  and  lag  flexures  so  that  the  first  flap  and  lag 
modes  essentially  involve  only  rigid-body  blade  motions.  Further,  the  design  ensures  a 
rotating  torsional  frequency  of  at  least  9/rev  over  the  entire  rotor  speed  range  tested  and, 
thus,  virtually  eliminates  the  need  to  consider  a  torsional  degree  of  freedom.  While  this 
model  is  not  representative  of  an  operational  rotor  system,  it  is  suitable  due  to  the  desire  to 
restrict  the  scope  of  the  physical  phenomena  being  studied  for  imtiai  research 
investigations  Future  experiments  should  introduce  the  blade  torsion  degree  of  freedom. 
The  experimental  damping  data  used  in  the  present  correlation  were  obtained  from  two 
rotor  models  of  identical  design,  Rotor  I  and  Rotor  II.  For  example,  at  the  collective  pitch 
angle  of  zero  degree,  the  damping  data  were  obtained  from  Rotors  I  and  II,  whereas,  at  the 
collective  pitch  of  three  and  six  degrees,  the  data  were  acquired  from  tests  on  Rotor  II. 
Procedures  were  essentially  the  same  for  both  the  rotors  Owing  to  manufacturing 
differences,  the  lag  and  torsion  stiffnesses,  nonrotating  flap  and  lag  frequencies  and 
structural  flap  damping  of  the  rotors  differ  slightly,  but  these  differences  are  too  small  to 
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be  significant.  Thus,  for  the  analysis,  only  the  Rotor  I  properties  were  used  (Table  1), 
since,  the  bulk  of  higu  advance  ratio  (p  >  0.3)  and  high  thrust  (CT/as  >  0  1)  data  refers  to 
model  Rotor  !(dotted  lines  in  Fig.  1). 

The  model  tested  was  a  63.84-in  diameter  three-bladed  rotor  mounted  on  a  very  stiff 
rotor  stand  so  that  the  stability  data  were  representative  of  an  isolated  rotor  The  rotor 
had  no  cyclic  pitch  control,  and  collective  pitch  was  set  manually  prior  to  each  run  At  a 
given  advance  ratio,  rotor  speed  and  collective  pitch,  the  shaft  tilt  was  the  only  means  of 
controlling  the  rotor.  Thus,  the  rotor  was  operated  untrimmed  with  a  virtually  unrestricted 
tilt  of  the  tip-path  plane.  The  lack  of  cyclic  blade  pitch  control  prevents  trimming  the  rotor 
to  arbitrarily  chosen  thrust,  propulsive  force,  Dapping,  or  hub  moment  conditions,  but  it 
does  remove  a  potentially  serious  source  of  friction  damping  from  the  system  and  ensures 
high  quality  data  (Ref.  1).  The  model  was  shaken  in  roll  to  excite  the  lag  modes  and  then 
locked  up  to  record  the  transient  response.  The  frequency  and  damping  data  were  then 
obtained  from  the  time  histories  via  a  moving  block  technique.  In  forward  flight  two  rotor 
speeds  were  used  (fl  =  750  and  1000  rpm),  corresponding  to  u,  values  of  0  72  and  0  6i 
Advance  ratio,  shaft  tilt  and  collective  pitch  were  varied  to  cover  the  test  envelope  shown 
in  Fig  1.  At  each  condition  tested,  at  least  two  separate  damping  measurements  were 
obtained. 


3.  Analytical  Model 

The  mathematical  model  represents  the  structural  simplicity  of  the  experimental 
model  rotor,  as  an  offset-hinged,  ngid  blade  executing  lag(inplane)  and  flap  motions  with 
coincident  spring  restrained  flap  and  lag  hinges.  The  hinge  offset  of  11  1%  is  accounted  for, 
in  both  the  rotor  steady  state  and  stability  analyses.  The  rotor  is  untrimmed,  with  no 
cyclic  pitch  (cyclic  flapping  present).  At  a  given  advance  ratio,  the  collective  pitch  and 
shaft  tilt  angles  are  the  known  control  input  parameters.  The  principal  ciastic  axes  of  the 
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hinges  do  not  rotate  with  biade  collective  pitch  angle. 

3.1  Dynamic  Stall  Modeling 

The  aerodynamic  representation  is  based  on  the  ONERA  dynamic  stall  models  of 
unsteady  lift  and  drag,  which  are  virtually  independent  (Refs  6,7).  With  T  representing 
the  instantaneous  circulation  providing  lift,  and,  rd  representing  the  analogous  circulation- 
like  drag  counterpart,  the  stall  dynamics  begin  with  the  following  representation 

r  =  r.-r  r3  and  rd  =  rd,  +  rdj  (i) 

where  subscripts  1  and  2  indicate  the  linear  and  the  nonlinear  (stalled)  components, 
respectively.  For  the  local  coordinates  of  Fig.  2,  the  corresponding  dynamic  stall  lift 
components  are  modeled  according  to  Eqs  (2)  : 


Ly  =  u,  (r,-rrj)  +  Uy  (rd,+rd3)  +  l0  (2a) 

U  =  -uy  (r,+rj)  +  ux  (rd,+rdj)  (2b) 

L0  =  b  s  Uy  (2c) 

k  f  i  4-  X  Pi  =  A  (aUyCOsa  -f  UCi,)  +  S  b  i  (2d) 

#  r,  *  2d-xk  r,  +  ts’O  +  d7)  r2  =  -  ^(l-i-d5)  (  UAC, 

+  e  i'(UxC0SQ  -r  UySinQ)  ACZ 

■f  ek  (Uycoso  -  Uxsino)  ]  (2e) 


where  F.  and  f:  represent  circulation  per  unit  length  for  the  linear  and  stalled 
contributions  Moreover,  the  parameters  b,  c.  and  s  denote  the  airfoil  semi-churd,  l.near  lift 
curve  slope,  and  apparent  mass  term,  respectively  In  Eq.  (2d),  the  term  UC.c  accounts  for 
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nonzero  lift  at  zero  angle  of  attack  and  the  difference  in  maximum  lift  at  positive  and 
negative  angles  of  attack  of  the  XACA  23012  airfoil  section.  The  term  <  signifies  the  airfoil 
rotation  rate  with  respect  to  air  mass  Therefore,  i  must  include  all  geometric  rotations  of 
the  airfoil,  even  if  they  come  from  sources  other  than  blade  pitch.  For  the  rigid  flap-lag 
motion  of  a  rotor  blade,  t  =  6  +  (1  +  <)  sin#  In  the  absence  of  cyclic  pitch  0  is  zero, 
hence,  c  =  (i  +  <)  sin5.  Further,  in  Eqs.  (2d)  and  (2e),  k  reflects  an  average  of  the  reduced 
frequency  of  the  free-stream,  and  coefficients  A,  6,  ts,  d,  and  e  are  the  stall  system 
parameters.  As  shown  in  Fig.  3,  ACt  is  the  difference  between  the  extrapolated  linear  lift 
coefficient  and  quasisteadv  stall  lift  coefficient  in  the  stalled  region  The  quasisteadv  stall 
lift  coefficient,  CM,  is  approximated  by  the  following  functions  for  angles  of  attack  between 
0°  and  180° 


Cis  —  C*i 

0°  <  a  <  a„° 

(20 

Cjs  =  sinQjj  cosqsj  +  CiQ 

o»c  <  a  <  45° 

(2g) 

Cts  ~  a  sina^j  cosa**  sin2a  +  Ci0 

45°  <  a  c  135° 

(2h) 

C zs  ~  —  a  sina**  cosa5S  +  Cj0 

135°  <  Q<(lS0-oa)° 

(2i) 

C«  =  Csi 

(lSO-a^C  a<  180° 

(2j) 

In  Eqs.  (2e)  to  (2i),  aM  is  the  quasisteady  stall  angle,  and  C,j  the  extrapolated  linear  lift 
coefficient,  which,  following  the  unified  lift  theory,  is 

Cti  =  a  sina  cosa  +  Ci0  0°<a<lS0°  (2k) 

To  ascertain  the  lift  coefficient  values  at  negative  angles  of  attack,  a  phase  shift  of  180*’  can 
be  applied  to  a.  It  is  noted  that  the  simplified  quasisteady  stall  model  incorporates 
different  maximum  lift  values  at  positive  and  negative  stall  conditions  as  expected  fur  the 
cambered  XACA  23012  airfoil  section,  however,  it  does  not  include  different  positive  and 
negative  stall  angles  of  attack. 


As  seen  from  Eqs.  (2c)  to  (2e),  there  are  six  coefficients  of  the  dynamic  stall  lift 
model,  identified  from  wind-tunnel  measurements,  and  for  the  OA212  airfoil  they  are  well 
discussed  in  Ref.  6  It  is  instructive  that  £  is  the  only  parameter  m  the  dynamics  equation 
of  the  unstalled  circulation.  Eq.  (2d),  that  is  dependent  on  ACZ.  In  a  general  sense 
applicable  to  a  series  of  airfoils  including  N’ACA  23012,  £  can  be  approximated  as 

£=  oCu/da-O.oef  1  +  1.43  |AC*|  )  (21) 

In  the  above  equation  0.5a  term  gives  the  same  instantaneous  step  response  (normalized 
lift),  one-half  the  final  value,  as  the  Wagner  (pitch)  function.  As  emphasized  in  Ref.  6, 
despite  the  differences  in  airfoils,  such  responses  of  normalized  lift  do  exhibit  qualitative 
similarity.  A  further  approximation  of  £  that  appreciably  simplifies  both  the  trim  and 
perturbation  analyses  is  to  remove  dependency  on  AC-,  since  £  is  expected  to  primarily 
influence  the  unstalled  or  linear  circulation  Ti,  as  typified  by  Eq.  (2d)  That  is, 

S=aCzi/da  -(a/2)  (2m) 

Fig.  4  shows  the  effect  of  this  approximation  on  the  dynamic  stall  lift  coefficient  for 
the  XACA  23012  airfoil,  the  experimental  data  is  from  Ref.  11  at  1  -  10s  Reynolds  number. 
Thus,  Eq.  (2m)  provides  a  viable  approximation  that  is  consistent  with  the  classical  linear 
unsteady  lift  theory.  Also,  it  removes  the  A'-'  dependency  of  the  unstalled  circulation  T 1. 
In  Fig.  4,  the  effects  of  further  simplification  in  £,  which  neglects  the  second  term  in  Eq. 
(2m),  are  also  shown.  It  is  seen  that  this  approximation  is  also  a  reasonable  alternative  for 
obtaining  a  qualitative  trend  However,  it  instantaneously  gives  the  normalized  lift 
response  that  is  equal  to  the  final  value  of  the  Wagner  function  and,  thus,  violates  the  time 
delay  in  the  classical  linear  unsteady  theory  (Ref.  6).  Hereafter,  £  according  to  Eq  (2m)  is 
used  in  the  dynamic  stall  lift  theory  unless  noted  otherwise  (see  Fig.  10,  which  explores 
sensitivity  of  lag  damping  to  such  changes  in  £). 

As  seen  from  Fig.  2,  the  resultant  velocity  U  and  the  angle  of  attack  a  are  given  by 
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U  =  ]  Uj  +  U*  =  j  U|  +  Uj  :  q  =  arctan(LyUx) 

(3a) 

UT  =  eh  cos<  +  cos$(l-K)  x+  ^sin(tHC) 

(3b) 

Xjp  =  -  eh  sin{  sin0  •fi^+Atcos°+/i  sinfl  cos(tM-{) 

(3c) 

Ux  =  UT  CO30O  +  Up  sin0o 

(3d) 

Uv  =  U?  sin0o  -  Up  cos 0O 

(3c) 

The  dynamic  stall  drag  model  follows  the  same  analytical  structure  of  the  lift 
model,  except  that  the  linear  component  Pa,  in  Eq.  (1)  does  not  involve  an  additional 
state.  The  dynamic  stall  drag  suitable  for  use  in  conjunction  with  the  dynamic  stall  lift 
model  is  formulated  according  to  the  following  equations: 

Ta,  =  U  Cdo  (4a) 

fa,  +  V  fa,  +  d  fa,  =  -  [  rj  U  ACrf  +  ed*  Uy  1  (*) 

where  AC  a,  as  shown  in  Fig.  5a,  indicates  the  difference  between  the  constant  profile  drag 
coefficient  and  the  quasisteady  stall  drag  coefficient,  analogous  to  AC:  ir.  Fig  3.  The  above 
formulation  of  dynamic  drag  in  terms  of  ctrculation-liVe  drag  variables  attempts  to 
incorporate  the  unsteady  wind-speed  variations  ana  the  heaving  motion  of  the  airfoil, 
analogous  to  the  circulation  lift  model. 

For  substall  conditions  with  positive  angles  of  attack  (e  g  hover  data  in  Fig.  1), 
Refs.  1  and  12  suggest  the  following  expression  for  the  quasisteady  stall  drag. 

Cds  =  0  0079  *  1  7  q1  (4c) 


According  to  Ref.  1,  Eq  (4c)  is  based  on  sieady-bending-momeni  data  from  a  test  with  the 
blades  similar  to  those  used  in  the  present  lest.  However,  the  present  test  also  includes 
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stalled  cases  in  forward  flight,  both  with  positive  and  negative  angles  of  attack  The 
quasisteady  stall  drag  data  appropriate  to  the  database  of  Fig.  1  are  not  available,  and  the 
present  study  uses  the  following  equation  : 

Cds  =  1.03 - (1.03-Cdo)  cos2(a+ae),  -180°  <  o<  180°  (4d) 

As  seen  from  Fig.  5, there  are  quantitative  differences  between  Eqs.  (4d)  and  (4c)  for 
moderate  angles  of  attack.  But  over  the  large  angle  of  attack  of  range  Eq  (4d)  agrees  well 
with  Eq.  (4c),  behaves  smoothly  for  a  complete  sweep  of  angle  of  attack,  and  is,  thus, 
expected  to  be  adequate  in  predicting  qualitatively  the  principal  effects  of  quasisteadv  stall 
drag  on  lag  damping. 

For  the  ONERA  OA212  airfoil,  the  dynamic  stall  lift  parameters  are  well  discussed 
in  Ref.  6.  By  comparison,  the  dynamic  stall  drag  model  is  still  evolving,  and  its  parameters 
ir.  Eq.  (4b)  are  assumed  to  be: 

ad  =  0.32  ,  rd  =  0.2  +  0.1  AC| ,  and  ed  =  -0.015  AC l  (5) 

3.2  Blade  Equations 

To  implement  the  dynamic  stall  lift  and  drag  effects  according  to  Eqs  (2)  and  (4), 
me  blade  is  discretued  into  X  radial  aerodynamic  elements,  and  the  aerodynamic  forces  are 
assumed  to  be  constant  in  each  element.  Accordingly,  the  rigid  flap-lag  equations  take  the 
form  : 

fl  +  Cq0  +  sin5cos5(l-r():  +  (4sm/fcos(  +  (-p-l)5  + 

-3siv*» 

i  *  ( 


(6a) 
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ccs3/)  (  +  (  -  2$\n0cos0(\+'Q0  +  t4  cos/5sin(  +  IV(  +  ZB 

s 

=  h  ms0  XFCs  r‘ Ai  (6b) 

S«  1 

In  Eqs.  (6),  and  F^  are  the  dimensionless  airloads,  acting  at  the  center  of  the  ith 

dement,  perpendicular  to  the  blade  and,  respectivdv,  perpendicular  and  paralld  to  the 
plane  of  rotation.  Moreover,  i-,  represents  the  radial  location  of  the  center  of  the  ith 

dement  from  the  blade  hinge,  and  Aj  is  the  length  of  the  ith  dement.  The  stiffness 

parameters  P,  IF  and  Z  in  Eqs.  (6)  are  defined  in  Ref.  2.  For  this  study,  P  =  1+u^  ,  IF  = 
ui  ,  and  Z  =0. 

According  to  Eqs.  (2),  (3)  and  (4),  the  Dap  and  lag  aerodynamic  forces  for  a  generic 
ith  blade  dement  can  be  expressed  as 

Fa  =  [  (r.+Tj)  UT  -  Up  (I\j,+rd,)  +  L0  cosl?0 1.  (7a) 

F0  =  -  ( (r.t-r-,)  up  +  uT  (rd,+rdJ)  +  ia  si  n0a }.  (7b) 

3.3  Solution  Scheme 

Concerning  trim  and  Floquet  eigenanalvsis,  a  few  details  are  added  here  for 
completeness.  The  rotor  is  operated  with  specified  control  inputs  for  collective  pitch  angle 
0o  and  shaft  tilt  as.  Since  cyclic  pitch  inputs  were  not  used,  the  rotor  is  not  trimmed  to  a 
specified  force  or  moment  condition.  The  equilibrium  conditions  are  found  by  time 
integration  over  several  rotor  revolutions  until  convergence  or  response  periodicity  is 
achieved.  Determination  of  uniform  steady  inflow  and  thrust  is  based  on  the  momentum 
theory.  It  must  be  emphasized  that  the  test  conditions  deviate  significantly  from  typical 
hdicopter  rotor  operations  (e.g.  large  negative  thrust  and  coning,  shaft  tilt  as  high  as  20® 
and  advance  ratio  varying  from  zero  to  0.55).  The  stability  margins  or  damping  levels  are 
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computed  from  a  standard  eigenanalysis  of  the  Floquet  transition  matrix  (FTM),  which 
was  generated  by  a  single-pass  initial  value  code.  Computational  reliability  is  verified  on 
the  basis  of  eigenvalue  condition  numbers  and  residual  errors  of  eigenpairs  (Ref  13) 

A  consist,  nt  perturbation  scheme  requires  that  all  quantities  depending  on  angle  of 
attack  be  perturbed  about  a  periodic  forced  response  with  equilibrium  angle  of  attack  a 
For  example,  with  a  =  a  -  a,  the  corresponding  perturbations  in  the  circulation  of  the  itt> 
element,  A Cr  and  AC<j,  are 


£,;  =  Tn-f,,  (8a) 

r,i  =  r  jj  -  r  3i  (sb) 


Though  the  required  perturbation  per  se  is  fairly  routine,  Eqs.  (8a)  and  (8b)  are  mentioned 
to  emphasize  the  roie  of  individual  discretized  elements,  and  Eqs  (8c),  (8d)  and  (8e)  are 
mentioned  to  emphasize  the  role  of  quasisteady  stall  lift  and  drag  characteristics  It  alsc 
should  be  emphasized  that  the  dynamic  stall  lift  parameters  f,  is,  d  and  e  in  Eqs  (2),  and 
dynamic  stall  drag  parameters  and  sd  in  Eq.  (5)  depend  upon  the  angle  of  attack  and 
require  perturbation. 


4.  Results 


If  not  stated  otherwise  the  results  in  Figs.  6  to  17  are  based  on  she  baseline  values 
given  in  Table  2,  and  specific  additional  values  are  so  identified  in  respective  figures  The 
lag  damping  predictions  are  based  on  the  uniform  steady  inflow  The  following  four  theones 
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are  examined.  1)  h near  theory,  a  constant  airfoil  lift-curve  slope  and  profile  drag  coefficient 
with  no  reverse  flow  (Figs.  6-16)  and  with  reverse  flow  (Fig  17),  2)  quasisteady  stall 
theory,  linear  theory  refined  to  include  quasisteady  nonlinear  airfoil  lift  and  drag 
characteristics  which  include  reverse  flow  effects  (Figs.  3  and  5),  3)  dynamic  stall  lift 
theory,  quasisteady  stall  theory  refined  to  include  dynamic  stall  lift,  and  4)  dynamic  stall 
theory,  which  includes  dynamic  stall  lift  and  dynamic  stall  drag  characteristics  of  the 
airfoil.  Given  the  extensive  scope  of  the  database  represented  in  Fig.  1,  only  a  portion  of 
the  correlation  is  presented. 

4.1  Forced  Response 

Figs.  6  and  7  show,  respectively,  the  trim  results  of  periodic  Gap-response  3  and 
angle  of  attack  a  over  one  rotor  revolution.  All  three  collective  pitch  angles  of  the  database 
in  forward  flight  (also  see  Fig.  1)  are  considered  for  an  advance  ratio  of  0  35  and  shaft  tilt 
angle  of  20°.  While  Fig.  6  shows  the  flapping  response  from  the  linear,  quasisteady  stall 
and  dynamic  stall  theories,  Fig.  7  shows  the  angle  of  attack  results  from  only  the  dynamic 
stall  theory  for  ihe  three  outer  blade  dements.  For  brevity,  ct  predictions  from  the  linear 
and  quasisteady  stall  theories  for  the  entire  blade  are  omitted  Similarly,  d  predictions 
from  the  dynamic  stall  theory  for  the  two  inner  segments  are  omitted  as  well,  although  the 
inner  segments  experience  much  higher  angles  of  attack,  they  are  aerodynamicallv  less 
effective.  As  seen  from  Fig.  6,  both  dynamic  stall  and  quasisteady  stall  effects  increase  with 
decreasing  collective  pitch  for  this  high  shaft  tilt  case  with  o»  =  20°.  Also,  as  stall  effects 
increase  the  linear  theory  p  results  deviate  significantly  from  those  of  the  stall  theories 
The  a  predictions  in  Fig.  7  support  the  3  predictions  in  Fig.  6,  in  which  the  maximum 
magnitude  of  the  local  angle  of  attack  a  decreases  substantially  as  collective  pitch 
increases.  Another  feature  merits  special  attention  .  most  of  the  blade  experiences  negative 
angles  of  attack  a  with  particularly  large  magnitudes  for  the  zero  collective-pitch  case.  This 
is  a  result  of  the  test  conditions  chosen  for  the  experimental  modd,  which  at  zero  collective 
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pitch  angle  and  at  high  advance  ratios  experienced  large  negative  angles  of  attack  and 
coning  angles  (Ref.  1)-  Since  the  blade  airfoil  is  cambered,  these  negative  angles  of  attack 
stall  conditions  represent  an  unusual  operating  condition  for  the  airfoil  For  the  higher 
collective  pitch  angles  of  3°  and  6°,  the  stall  regions  shrink  noticeably  on  the  retreating 
blade,  particularly  near  the  blade-tip  where  a  essentially  remains  unstailed,  e  g  see  Figs 
7b  and  7c  for  the  outermost  fourth  and  fifth  elements. 

The  dynamic  stall  effects  on  the  trim  results  are  further  pursued  in  Fig  8,  in  which 
thrust  levels  C T/<rs  based  on  the  linear  theory  are  shown.  Results  for  dynamic  stall  theory 
(not  shown)  indicated  that  the  linear  theory  overpredicted  C T/<rs  by  about  25%  for  0o  =  C° 
at  it  =  0  35.  Nevertheless,  Fig.  8  serves  as  a  rough  guideline  in  illustrating  the 
interdependence  of  thrust  and  stalled  conditions  For  example,  comparison  of  Figs  C  and  8 
shows  that  the  quasisteady  stall  and  dynamic  stall  effects  arc  more  pronounced  at  the  high 
thrust  conditions.  Regarding  the  database  in  Fig.  1,  the  tnm  results  in  Figs  6,  7  and  8 
show  that  dynamic  stall  should  play  a  major  role  for  the  0°  collective  case,  particularly  for 
the  high-advance-ratiu  and  high-shaft-tilt  combinations.  That  role  should  be  minor  for 
0o  =  3'  and  6°,  the  higher  the  collective  pitch,  the  less  the  dynamic  stall  effects 

4.2  Lag  Damping  Correlation 

The  following  damping  predictions  are  based  on  the  nonlinear  quasisteady  drag 
characteristics  using  Eq.  (4d).  Fig.  S  shows  lag  damping  versus  rotational  speed  in  hover 
for  collective  pitch  angles  60  =  4s  and  8°.  Here,  the  dynamic  stall  effects  are  negligible,  as 
determined  from  the  stall  angles  based  on  the  tnm  values  ([angle  of  attack  aj  <  12°)  The 
predictions  cluster  into  two  groups  of  linear  theory  vis-a-vis  the  other  three  thcones  The 
linear  theory  is  fairly  adequate  in  the  prediction  of  damping  level  trends,  including  the 
reduction  in  damping  level  around  375  rpm,  where  the  flap  and  lag  frequencies  coalesce. 
However,  the  discrepancies  between  the  data  and  the  linear  theory  consistently  increase  for 
n  >  450  rpm.  By  comparison,  the  remaining  three  theones  provide  excellent  correlation 
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throughout.  The  difference  between  the  three  stall  theories  and  the  linear  theory  is 
primarily  due  to  the  nonlinear  drag  effects.  Furthermore,  the  difference  between  the 
quasisteady  stall  theory  and  the  dynamic  stall  theories  is  rather  small  overall,  but 
consistently  increases  with  increasing  collective  pitch  and  rotational  speed,  it  illustrates  the 
the  linear  unsteady  lift  effects,  since,  stall  is  r.ot  an  issue  here. 

The  lag  damping  correlation  in  forward  flight  is  treated  in  Figs.  10  to  17  for  a 
complete  range  of  advance  ratio  (0  <  /i<  0.55)  and  collective  pitch  angle  (0o  —  0°,  3°  and 
6°),  with  the  shaft  tilt  angle  Qs  varying  from  0°  to  20°.  Fig.  10  shows  the  effects  of 
simplifying  6  according  to  Eq  (2m)  on  lag  mode  damping,  and,  thus,  it  complements  the 
results  in  Fig.  4.  The  damping  is  presented  versus  advance  ratio,  the  shaft  angle  as  =  16° 
for  both  collective  pitch  angles  of  0°  and  3°,  and  as  =  20“  for  collective  pitch  angle  of  6° 
The  predictions  are  based  on  the  dynamic  stall  lift  theory  for  two  representations  of  £ 
according  to  Eq.  (21),  and  simplification  according  to  Eq.  (2m).  With  increasing  advance 
ratio,  both  predictions  show  the  same  trend  of  an  initial  reduction  in  the  damping  followed 
by  an  abrupt  increase  in  the  damping.  That  this  trend  is  not  sensitive  to  such  changes  in  £ 
is  an  appealing  feature  of  the  dynamic  stall  lift  model. 

Fig.  11  shows  lag  damping  versus  advance  ratio  for  shaft  tilt  a*  =  0“  and  4°,  both  at 
0“  collective  pitch  angle.  For  these  shaft  tilt  angles,  the  thrust  levels  ICT/o»|  are  very  low 
for  all  advance  ratios  with  a  maximum  value  of  about  0.03  (see  Fig  8a)  This  is  well 
reflected  by  the  almost  identical  predictions  from  all  four  theories  for  M  <  0  4  At  advance 
ratios  above  0.4,  the  linear  theory  predicts  higher  values  than  the  stall  theories,  and  this 
deviation  from  the  stall  theories  and  the  data  increases  with  increasing  advance  ratio  On 
the  other  hand,  the  quasisteady  stall  and  dynamic  stall  lift  theories  continue  t_  predict 
almost  identical  values  for  ^  >  0  4.  This  shows  negligible  effects  of  dynamic  stall  lift, 
including  unsteady  lift  in  substall  for  0o  =  0°.  Further  study  shows  that  the  deviation  of 
the  linear  theory  for  *i  >  0.40  is  associated  with  very  high  lift  prediction  due  to  the 
assumption  of  constant  lift  curve  slope  and  no  reverse  flow. 
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As  seen  from  Fig.  12,  the  lag  damping  correlation  is  continued  for  3°  collective  p'tch 
angle  for  three  relatively  low  values  of  shaft  tilt,  Qj  =  0°,  4°  and  8°.  Here  also  the  thrust 
level  |CTMj  is  relatively  low  and  it  generally  decreases  with  increasing  shaft  tilt  (see  Fig 
8b).  The  data  are  available  for  the  low  advance  ratio  range,  n  <  0.3,  and  all  four  theories 
correlate  well  with  data.  Thus,  for  n  <  0.3,  the  quasisteady  stall,  dynamic  stall  lift  and 
dynamic  stall  theories  predict  essentially  the  same  damping,  and  the  linear  theory  deviates 
from  these  three  theories  mainly  due  to  nonlinear  drag  effects  in  substall.  Thus,  at  low 
shaft  tilt  angles  treated  in  Figs  11  and  12,  the  quasisteady  stall  theory  and  the  two 
dynamic  stall  theories  provide  excellent  correlation. 

The  damping  correlation  at  relatively  high  shaft  tilt  angles,  a$  =  8°,  12°.  16°  and 
20°,  at  0o  =  0°,  is  given  in  Fig.  13.  Up  to  advance  ratio  of  about  0.2  (|CT/os|  <  0  06),  the 
predicted  damping  essentially  remains  constant,  all  four  theories  are  virtually  identical  and 
correlate  wdl  with  the  data.  For  advance  ratio  n  >  0.25,  the  data  show  that  the  damping 
increases  rather  abruptly,  this  abruptness  bang  more  marked  with  increasing  shaft  tilt  os 
and  advance  ratio  u.  Fig.  Sa  facilitates  a  better  understanding  of  this  abrupt  and  marked 
increase  in  damping.  It  shows  that  the  thrust  level  vanes  from  -0.07  at  qs  =  8°  to  -0.23  at 
Of  =  20°,  basically  representing  high  negative  thrust  and  highly  stalled  conditions. 
Surpns-.ngly,  the  linear  theory  predicts  this  trend  of  abruptly  and  markedly  increasing 
damping.  The  quasisteady  stall  theory  eventually  predicts  this  peculianty,  however,  the 
structure  or  pattern  of  this  increase  is  different  from  that  of  the  data,  especially  at  shaft 
tilt  angles  of  16°  and  20°  for  which  quasisteady  stall  theory  predicts  mild  instability  and 
then  shows  the  abrupt  increase.  But  inclusion  of  dynamic  stall  lift  effects  bnngs  in 
significant  quantitative  improvement  in  the  correlation,  more  or  less  qualitative 
improvement  at  high  shaft  tilts  of  IS0  and  20°.  Inclusion  of  dynamic  stall  drag  yields 
further  quantitative  improvement.  It  must  be  mentioned  that  these  improvements  of 
dynamic  stall  lift  and  iynamic  stall  drag  are  consistent  in  that  the  .mprovements  increase 
markedly  with  increasing  advance  ratio  It  must  be  noted  however,  that  overall,  the  linear 
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theory  gives  the  best  correlation  for  the  case  of  zero  collective  pitch  angle. 

The  correlation  for  3°  collective  pitch  angle  is  given  in  Fig.  14  for  high  shaft  tilt 
angles,  as  >  12°.  The  data  show  the  clear  trend  of  slowly  decreasing  damping  with 
increasing  advance  ratio  However,  the  linear  theory  is  qualitatively  inaccurate  here, 
predicting  the  opposite  trend  of  rather  rapidly  increasing  damping  with  increasing  advance 
ratio.  The  other  three  theories  cluster  together  and  predict  the  trend  of  the  data  better.  A 
significant  observation  is  that  the  three  stall  theories  for  high  p  and  a5  combinations 
(sav  p  >  0  4  and  as  >  14)  predict  an  abruptly  increasing  damping,  as  do  the  linear  theory 
in  general  and  the  data  in  Fig.  13  for  0a  =  0C  in  particular.  In  the  absence  of  experimental 
data  at  high  p  for  0o  =  3°,  it  is  difficult  to  judge  the  validity  of  the  dynamic  stall  theory 
here. 

The  correlation  for  the  6°  collective  pitch  angle  is  shown  in  Fig.  15  for  os  =  20°  and 
is  quite  similar  to  that  in  Fig.  14  where  the  damping  slowly  decreases  with  increasing 
advance  ratio.  This  trend  is  weil  predicted  by  the  quasisteady  stall,  dynamic  stall  lift  and 
dynamic  stall  theories.  In  contrast,  the  linear  theory  is  qualitatively  inaccurate.  The 
improvement  over  the  linear  theory  by  the  quasisteady  stall  theory  is  significant  and 
reflects  the  nonlinear  drag  effects  on  lag  damping.  The  dynamic  stall  theories  deviate 
marginally  from  the  quasisteady  stall  theory  mainly  due  to  the  unsteady  linear  lift  effects. 
With  increasing  advance  ratio,  particularly  for  p  >  0.4,  the  dynamic  stall  lift  and  dynamic 
stall  drag  come  into  play.  As  in  the  previous  case  for  0o  =  3°  in  Fig  14,  here  also  it  is  not 
possible  to  validate  such  high-advance-ratio  predictions  on  the  basis  of  test  data. 

To  facilitate  an  overall  assessment  of  the  adequacy  of  the  four  theories  for  a  wide 
range  of  data  in  forward  flight,  Fig.  16  shows  a  summary  of  the  correlation  for  the 
collective  pitch  angles  of  0°,  3°  and  6°  at  high  shaft  tilt  angles  of  14°,  14°  and  20°, 
respectively.  The  predictions  of  l.near  and  quasisteady  stall  theories  are  unsatisfactory, 
each  providing  good  correlation  for  one  or  two  cases  and  qualitatively  inaccurate  prediction 
for  the  remaining  cases.  For  example,  the  linear  theory  gives  good  correlation  for  the  case 
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of  zero  collective  pitch  angle.  Fig  16a.  But  for  3°  and  6°  collective  pitch  angles,  Figs.  16b 
and  16c,  it  predicts  increasing  damping  with  increasing  advance  ratio  which  is  qualitatively 
inaccurate  in  comparison  with  the  test  data  In  contrast,  quasisteady  stall  theory  gives 
good  correlation  for  collective  pitch  angles  0o  =  3°  and  6°.  However,  for  0o  =  0°, 
quasisteady  stall  theory  is  both  quantitatively  and  qualitatively  inaccurate.  On  the  other 
hand,  the  dynamic  stall  lift  theory  provides  the  crucial  correction  for  0o  =  0“  and  gives 
fairly  accurate  prediction  in  all  three  cases.  In  other  words,  quasisteady  stall  theory  refined 
to  include  dynamic  stall  consistently  provides  the  best  correlation  for  the  broadest  range  of 
experimental  data.  Furthermore,  it  should  be  noted  that  the  linear  theory  in  highly- 
separated  flow  conditions,  for  the  collective  pitch  angle  of  zero  degree  at  high  shaft  tilts, 
provides  surprisingly  adequate  correlation. 

The  summary  results  in  Fig.  16  may  also  be  interpreted  in  terms  of  the  variation  of 
thrust  level  with  advance  ratio.  In  Fig.  16a,  C T/os  varies  from  roughly  0  to  -0.1  over  the 
advance  ratio  range  of  0  to  0.4.  In  Fig.  16b,  CT/cr$  varies  from  0.025  to  -0.045,  and  in  Fig. 
16q  from  0.05  to  -0.05  Therefore,  Fig  16a  changes  from  low  thrust  to  high  negative  thrust 
with  significant  airfoil  stall.  Fig.  16b  changes  from  moderate  positive  thrust  to  moderate 
negative  thrust  with  moderate  stall,  and  Fig.  16c  changes  from  high  positive  to  moderate 
negative  thrust  with  similarly  moderate  airfoil  stall. 

With  respect  to  the  effect  of  rotor  thrust  level  and  advance  ratio  on  theoretical 
predictions  of  lead-lag  damping,  the  following  observations  can  be  made.  For  both  linear 
and  nonlinear  theories,  the  basic  trend  is  for  damping  to  increase  with  advance  ratio  after 
an  intermediate  reduction  in  damping  from  the  nonlinear  theories.  This  basic  trend  can  be 
interpreted  as  the  result  of  two  effects  known  to  influence  the  lead-lag  damping.  1)  the 
effect  of  thrust  variation  with  advance  ratio  and  2)  periodic  coefficient  effects  due  to 
advance  ratio  The  basic  difference  between  the  linear  and  nonlinear  theory  is  that  the 
nonlinear  theory  predicts  that  the  damping  decreases  before  the  increased  trend  takes 
place.  This  initial  decrease  in  the  damping  is  related  to  the  reduced  lift  prediction  by  the 
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stail  theories.  Comparing  the  nonlinear  stall  theories,  it  is  interesting  that  the  effect  of 
dynamic  stall  in  forward  flight  tends  to  offset  the  large  influence  of  quasistatic  stall  on 
lead-lag  damping,  although  this  compensation  is  small  for  hover  and  low  advance  ratios. 

From  Fig.  16a  (also  Figs.  13),  it  is  seen  that  the  earlier-used  conventional  linear 
theory,  Ci  x  a  a,  provides  a  very  satisfactory'  correlation  at  0o  =  0°.  This  adequacy  of  the 
linear  theory  for  the  highly  stalled  conditions  is  surprising  and  merits  reassessment.  As 
noted  earlier,  this  linear  theory  assumes  a  constant  iift-curve  slope  which  predicts 
unrealistically  high  lift  at  large  incidence  and,  also  excludes  reverse  flow  effects,  an 
exclusion  that  further  augments  over-prediction  of  lift.  The  linear  theory  was  refined 
according  to  Eq.  (2k),  Ci  *  a  sino  costs,  to  account  for  the  large  angles  of  attack  and 
reverse  flow  effects,  henceforth  referred  to  as  the  refined  linear  theory.  Accounting  for 
reverse  flow  according  to  Eq.  (2k)  may  not  be  the  most  accurate  way,  but,  it  does 
incorporate  the  correct  sign  reversal  in  the  lift  and  also  predicts  a  more  realistic  lift 
(see  Fig.  3).  Fig.  17  compares  the  the  refined  linear  theory  with  the  linear  theory,  which 
assumes  a  constant  lift-curve  slope  and  no  reverse  flow.  While  the  solid  lines  in  Figs.  17 
correspond  to  the  refined  linear  theory,  the  broken  lines  correspond  to  the  linear  theory 
without  reverse  flow  effects,  as  used  in  Figs  9  to  16.  The  three  cases  are  identical  to  those 
in  Fig.  16.  It  is  seen  from  Fig  17  that  the  refined  linear  theory  qualitatively  alters  the 
damping  predictions.  For  example,  at  zero  collective  pitch,  Fig  17a  shows  that  the  refined 
linear  theory  fails  to  capture  the  increasing  damping-level  trend  of  the  test  data  for  high 
advance  ratios.  In  fact,  it  shows  decreasing  damping  in  contrast  to  the  predictions  from  the 
linear  theory  without  reverse  flow.  Fig.  17b  also  shows  the  significant  effects  of  large 
incidences,  but,  significantly  enough,  the  refined  linear  theory  predictions  are  closer  to  the 
test  data.  For  the  case  of  0o  =  6°  at  shaft  tilt  of  20°,  Fig.  I7c  shows  that  both  the  versions 
of  the  linear  theory  fail  to  capture  the  trend  of  the  test  data  adequately.  Overall,  Fig.  17 
demonstrates  that  even  the  refined  linear  theory  is  not  viable  at  high  advance  ratios  when 
the  blade  encounters  high  angles  of  incidence.  Moreover,  rhe  seeming  adequacy  of  the  linear 
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theory  without  reverse  flow  for  the  isolated  case  of  zero  degree  collective  pitch  is  erroneous. 

The  relative  prediction  capabilities  of  the  various  theories  may  be  summarized.  The 
linear  theory  is  successful  when  the  rotor  thrust  level  is  low,  i.e.,  when  both  the  collective 
pitch  and  advance  ratio  are  low  (as  in  Fig.  16a)  or  when  a  combination  of  high  collective 
pitch  offsets  a  high  shaft  angle  (as  in  Fig.  16c  at  an  advance  ratio  near  0  25).  The 
quasisteady  stall  theory  is  adequate  when  thrust  level  is  moderately  high  and  advance  ratio 
is  low  (as  in  hover,  Fig.  9,  and  Fig  16c  at  low  advance  ratio).  It  is  interesting  to  note  that 
a  continuous  shift  between  these  two  conditions  is  exhibited  m  Fig.  16c  where  quasisteady 
stall  successfully  predicts  the  increased  damping  compared  to  linear  theory  for  high  thrust 
m  hover,  but  the  two  predictions  begin  to  merge  near  an  advance  ratio  of  0.15  when  airfoil 
stall  effects  vanish  as  the  thrust  level  drops  to  zero.  Finally,  the  dynamic  stall  theory  is  the 
most  effective  theory  when  both  the  thrust  level  and  the  advance  ratio  are  large. 

5.  Conclusions 

The  detailed  analytical  investigations  and  correlation  of  lag  damping  with  the 
experimental  data  lead  to  the  following  major  findings  : 

1)  Conventional  Floqucl  analysis  with  ONERA  unified  dynamic  stall  lift  and  drag  models 
provide  a  viable  approach  for  predicting  lag-mode  damping  in  stalled  forward  flight 
conditions. 

2)  Tor  low-advance-ratio  cases  (p  <  0  20),  the  nonlinear  substaj]  lift  and  mainly  drag 
characteristics  of  the  airfoil  section  significantly  improve  the  correlation  with  the  low- 
Reynolds-numter  model  lest  data,  confirming  earlier  analytical  and  experimental 
investigations  restricted  to  hovering  conditions 

3)  The  linear  theory  with  constant  lift-curve  slope  is  not  reliable  except  when  the  rotor 
thrust  is  low  Its  seeming  adequacy  for  the  collective  pitch  of  zero  degree  is  erroneous,  due 
to  neglecting  the  effects  of  large  angles  of  attack  and  reverse  flow  Overall,  the  refined 


linear  theory,  accounting  for  the  large  angles  of  attack  and  reverse  flow  effects, 
qualitatively  differs  from  the  linear  theory  without  refinements,  but  it  is  also  found  to  be 
inadequate  except  when  the  rotor  thrust  is  low. 

4)  The  quasisteadv  stall  theory  is  a  significant  improvement  over  the  linear  theories, 
without  and  with  refinements,  when  the  advance  ratio  is  low.  However  it  overpredicts  the 
effect  of  stall  at  higher  advance  ratios. 

5}  Overall,  the  dynamic  stall  lift  theory  provides  significant  and  consistent  quantitative 
and  qualitative  improvement  over  the  quasisteady  stall  theory  and  is  found  to  be  viable  for 
the  entire  database 

6)  Dynamic  stall  drag  generally  provides  further  quantitative  improvement. 

7)  The  lack  of  sufficient  experimental  data  for  high  thrusi,  and  moderateto-high  advance 
ratio  conditions  limits  the  ability  to  fully  evaluate  the  complex  character  of  the  prediction 
methods.  Furthermore,  the  low  Reynolds  number,  and  the  negatively  stalled  conditions  of 
the  cambered  airfoil  should  be  removed  in  future  experimental  investigations. 
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Table  1  :  Experimental  Model  Rotor  I  Properties 


Number  of  blades  3 

Airfoil  section  N'ACA  23012 

Hover  blade-tip  Mach  number  at  1000  rpm  0.25 

Hover  blade-tip  Reynolds  number  at  1000  rpm  240,000 

Rotor  radius,  m  31.92 

Blade  chord,  in  1.65 

Radial  location  of  the  center  of  flexures,  in  3.55 

Radial  location  of  the  blade  center  of  mass,  in  11.10 

Mass  outboard  of  the  flexure  center,  slugs  0.013 

Flap/lag  inertia  about  flexure  center,  sfogs-tn 1  1.8 

Torsional  inertia,  sfeys-in’  0.00294 

rlexure  assembly  flap  stiffness,  in-lb/ red  66.0 

Flexure  assembly  lag  stiffness,  in-fl/  rod  281.0 

Flexure  assembly  torsional  stiffness,  in- lb/ rod  215.0 

Nomotating  flap  frequency,  Hz  3.09 

Nonrotating  lag  frequency,  Hz  6.98 

Nonrotating  torsion  frequency,  Hz  149 

Nonrotating  flap  damping,  stc'1  -  0.15 

Nonrotating  lag  damping,  sec'1  -  0  09 

Blade  pretwist 


0.0 
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Table  2  :  Baseline  Parameters 


Lift  curve  slope,  a,  rrd'1 

Lock  number,  7  (based  on  a=  5.73  ) 

ProGle  drag  coefficient,  C<j0 

Camber  of  NACA  23012  airfoil,  Oc 

Lift  coefficient  at  <r=0c;  Ci0 

Quasisteadv  stall  angle  of  NACA  23012  airfoil,  a* 

Rotor  solidity,  as 

Nondimensional  hinge  offset,  eh 

Nondimensional  blade  root  cutout 

Rotor  speed,  fl 

Number  of  blade  segments,  N 

Apparent  mass  quantity,  s 


5.73 
7.54 
0.0079 
1.5° 
0.15 
±12° 
0.0494 
0.111 
0.186 
1000  rpm 
5 

0.0 


Other  structural  parameters  corresponding  to  the  experimental  rotor  are  given  in  Table  1. 
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F’igure  2  :  Sectional  aerodynamic  forces  and  velocities 


Fig.  3  :  Quasisteady  and  Linear  Lift  Coefficient  versus  Angle  of  AAtack 
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Fig.  4  :  Effect  of  Stall  Parameter  6  on  Lift  Hysteresis 
Loops  of  NACA  23012  Airfoil 
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Fig.  5  :  Quasisteady  and  Linear  Drag  Profiles 
for  Low  Reynolds  Number  Flow 
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Fig.  S  Periodic  Flap  Response  in  Forward  Flight.  =  0.35.  a ,  =  20° 
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Fig.  7  Equilibrium  Local  Angles  of  Attack  Variations,  n  -  0  35.  a,  =  2i 
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Fig.  1 1  :  Lag  Damping  Correlation  in  Forward  Flight 
for  Low  Shaft  Angles,  &a  =  0° 
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Fig.  12  :  lag  Damping  Correlation  in  Forward  Flight 
for  low  Shaft  Angles.  d0  -  3° 
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Fig.  14  •  Lag  Damping  CorrtlaUon  in  Forward  iTight 
for  High  Shaft  Angles,  -  3s 
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Fig.  15  :  Lag  Damping  Correlation  in  Forward  Flight,  d0  =  6°,  as  =  20° 


sec“‘  a, 3cc“*  a.  sec 


FI*  16  Sucunarj  of  Uf  Datnptn*  Correlation  in  Forward  Flifbl 
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Fig.  17 :  Adequacy  of  Linear  Theory  in  Forward  Flight 
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